-->[f,xopt,gradopt,work]=optim(purequad,x0vec,imp=9)
Norm of projected gradient lower than   0.7321859D-34.

 work  =
 
 
         column 1 to 6
 
    2.  - 5.551D-17    6.    1.355D-20    5.421D-20    2.711D-20  
 
         column  7 to 11
 
    3.253D-19  - 1.355D-20  - 5.421D-20  - 2.711D-20  - 3.253D-19  
 
         column 12 to 15
 
    6.019D-36    1.204D-35    1.204D-35    7.222D-35  
 gradopt  =
 
   1.0D-34 *
 
    0.1203706    0.7222237  
 xopt  =
 
   1.0D-34 *
 
    0.0601853    0.1203706  
 f  =
 
    4.709D-70  
 
-->diary(0)